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I""! Abstract 

X 

O Measuring the muon flux is important to the Sanford Underground Labo- 

ratory at Homestake, for which several low background experiments are being 
planned. The nearly-vertical cosmic ray muon flux was measured in three loca- 
tions at this laboratory: on the surface (1.149±0.017 x 10~^ cm~^ s~^ sr"-'^), at 
the 800-ft (0.712 km w.e.) level (2.67 ±0.06 x 10"^ cm^^ g-^ sr^i), and at the 
2000-ft (1.78 km w.e.) level (2.56 ± 0.25 x lO""^ cm^^ g-^ sr^i). These fluxes 
agree well with model predictions. 
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C ■ The Homestake Mine in Lead, South Dakota, USA was identified in 2007 



as the final candidate site for the Deep Underground Science and Engineering 
Laboratory (DUSEL). In advance of the federal funds to further develop the 
mine, the South Dakota Science and Technology Authority (SDSTA), which 
operates the Sanford Underground Laboratory, is offering an early science pro- 
gram mainly to characterize aspects of the site environment. It is located at 



5—1 44.35° N, 103.77° W, with a surface elevation of 1620 m above sea level. Ini- 

^ tially, the LUX (Large Underground Xenon) dark matter search T and the 

Majorana neutrinoless double beta decay experiment ^2 will be located there. 
Measurements of external backgrounds, including cosmic ray muons as well as 
gammas and neutrons, will be of paramount importance to the design of these 
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sensitive rare-event searches. 

The differential muon flux at the 4850 ft. (4.40 km water equivalent) level 
at Homestake was measured by Cherry et al [3]. However, more measurements 
are needed to characterize the muon flux as a function of depth. This paper 
describes new measurements of the cosmic ray muon flux at three locations: the 
surface (in a building that provided ~1 m w.e. of shielding), the 800 ft. (0.712 
km w.e.) level, and the 2000 ft. (1.78 km w.e.) level. We present a comparison 
between our reported flux and the model of [4] , and we find agreement between 
our experimental results and that model. 

Cosmic-ray muon flux as a function of depth has been studied by many 
experiments at underground facilities around the world. Measurements of the 
muon flux per unit solid angle as a function of slant depth from Castagnoli [S] , 
Barrett [6], Miyake [7], WIPP [8], Soudan [9], Kamioka [lO], Boulby [II], Gran 
Sasso [121 [13], Frejus [H] and Sudbury [IS] have been used to develop a model, 
which can be used to predict the muon flux per unit solid angle as a function of 
depth [1]. 

1. Methods 

The muon detector consists of four fast plastic scintillation counters (Saint- 
Gobain BC-408), each a square with a side length of 30.5±0.1 cm. As shown in 
Figure [T] the distance from the top to the bottom counter is 64.0 cm. Two of 
these counters are 0.5 cm thick, and the others are 1.0 cm thick. Each is coupled 
by a trapezoidal acrylic lightguide to a Photonis XP2020 photomultiplier with 
Photonis VD124K base. Waveforms from each detector are recorded by a 12-bit 
flash analog-to-digital conversion module with a sampling frequency of 170 MHz; 
it filters the data onboard with field programmable gate arrays and transmits 
digitized pulses through an Ethernet interface to a standard personal computer. 

The detector station also includes a 1.2 L liquid scintillation counter filled 
with Eljen Technologies EJ-301 or EJ-309 material. This counter has been used 
to study techniques for neutron counting in the underground environment; a 
future publication will quantify neutron backgrounds in the mine. 

The gamma ray flux at each of the sites is at the level of 1 cm^^ s~^; the 
counting rate for background gamma events in each of the plastic scintillator 
detectors is therefore ^1 kHz. Nearly all of the gamma flux is at energies of 
2.5 MeV or less, as was reported in pB]. A substantial coincidence requirement 
is needed to distinguish cosmic ray muon events from the gamma background. 
At the 2000 ft. depth, a two-fold coincidence analysis was found to be domi- 
nated by the gamma background, and at shallower depths they still represent 
a substantial correction. Consequently, all of our results are based on studies 
where we require that at least three of the four detectors record a pulse within a 
~70 ns time interval. Because of the higher counting rate and reduced sensitivity 
to several systematic uncertainties, we present the three- fold coincidence anal- 
ysis as our primary result, and a four-fold coincidence analysis (where all four 
detectors are required to fire simultaneously) as a partially independent check. 
In particular, the agreement between the three-fold and four-fold coincidence 
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Figure 1: Vortical positioning of elements of the muon detector. 



results demonstrates that secondary particles such as high-energy electrons pro- 
duced by muon interactions do not affect the counting rate substantially relative 
to the precision of this measurement. 

A geometric Monte Carlo calculation was used to determine the solid angle 
accepted by each of these analysis methods. It took into account only the size 
and position of each of the detector elements, assuming straight muon tracks. 
As shown in Figure [2] it was used to determine the acceptance probability as 
a function of polar angle P{d), which was then integrated to determine the 
accepted solid angle: 

t/2 



n 



In I P{e) si 



sin e d6 



When a four-fold coincidence is required, the accepted solid angle is 0.226 sr; 
it is 0.457 sr when only a three-fold coincidence is required. If we assume an 
incident muon distribution proportional to cos^ 9, 90% of the flux in the three- 
fold coincidence analysis would be within 25° of vertical, while 90% of the flux 
in the four-fold coincidence analysis would be within 19° of vertical. 

This Monte Carlo program was also used to study the systematic uncertainty 
arising from misalignment of the detector elements. A horizontal displacement 
of one detector element in this program by 2 cm, which is believed to represent 
the worst realistic possibility for the data collected on the 800 ft. level, changed 
the calculated solid angle by a maximum of 0.4% for the three-fold coincidence 
analysis and 1.1% for the four-fold coincidence analysis. The alignment on other 
levels is believed to have been substantially better, with a maximum possible 
displacement of 1 cm. 

The energy scales of the detectors were calibrated based on the observed 
pulse amplitude spectra for four-fold coincident events. All such events on 
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Figure 2: Acceptance probability as a function of polar angle 9 for the four-fold coincidence 
analysis (solid line) or the three- fold coincidence analysis (dotted line). 



the surface and at the 800 ft. level were presumed to be minimum- ionizing 
particles that would give a most probable energy deposition of 1.0 MeV in 
the thin detectors and 2.0 MeV in the thick detectors. These energies were 
computed from the scintillator density of 1.032 g/cm'^, assuming the minimum- 
ionizing dE/dx given by the Bethe-Bloch equation 17J. At the surface and the 
800 ft. depth, there was sufficient statistical power to allow the calibration to 
be determined in situ. At the 2000 ft. depth, a calibration from the surface was 
applied; this method would have been preferred at the 800 ft. depth as well, 
but detector repairs required recalibration to be completed underground. 

The digitization hardware thresholds were set as low as practical, corre- 
sponding to approximately 0.4 MeV for the thin detectors and 0.9 MeV for the 
thick detectors. Analysis thresholds were then established at 0.75 MeV in the 
thin counters and 1.5 MeV in the thick counters. These thresholds require a 
total energy deposition of at least 3 MeV for a three-fold coincidence, which is 
beyond the endpoint of the gamma spectrum, but still maintains an efficiency 
for muons that can be determined effectively. 

The efficiency corrections associated with these energy cuts were determined 
from the data. At the surface, we assumed that all four-fold coincident events 
were caused by minimum-ionizing particles; other particles from atmospheric 
showers would have been shielded effectively by the ^1 m w.e. provided by the 
building above the detector. Therefore, we computed the ratio of the number of 
events that passed the analysis cuts to the the total number of recorded four-fold 
coincidences where only the hardware thresholds were used. 

Having determined this efficiency on the surface, we then applied it to data 
collected at the 2000 ft. level. The uncertainty associated with this procedure 
includes a statistical component from the finite number of coincidences at the 



surface. However, another part of the uncertainty is related to the stabiUty of 
the detector gain, and therefore the energy scale, especially as the detector was 
being moved underground. We checked for these gain shifts using three standard 
beta/gamma calibration sources, ^^Na, ^°Co, and ^°Sr, collecting pulse height 
spectra before and after relocating the detector. Averaging the results from 
these sources, we found gain reductions of 0.8%, 2.0%, 4.7%, and 3.1% for the 
four detector elements. We corrected the energy scale in the analysis according 
to these results, and we treated the 1.0% uncertainty in each detector calibration 
as the systematic error for the efficiency in the result from the 2000 ft. level. 
Examination of the energy distribution in each detector showed that the single- 
detector rate would change by 0.4% for each percent change in the assumed 
energy scale, and a Monte Carlo simulation verified the naive predictions that 
such a change would give a modification of 0.4% to the three-fold coincidence 
analysis and 1.6% to the four-fold coincidence analysis. 

Because it was necessary to recalibrate the detector while it was at the 800 ft. 
level, a different procedure was applied there. We determined separately for each 
scintillator the fraction of four-fold coincidence events where only the hardware 
threshold was imposed that would be cut by applying the energy threshold in 
that scintillator. We incorporated these efficiencies into the geometric Monte 
Carlo simulation program and observed the change in the number of accepted 
three-fold and four-fold coincidence events in order to determine the overall ef- 
ficiency. The uncertainty associated with this procedure is taken to be purely 
statistical; any systematic drift in the energy scale should present a negligible 
effect, as the detector remained untouched underground at nearly constant tem- 
perature. The computed efficiency was substantially higher than before, because 
an upgraded flash ADC electronics module was installed when the detector was 
relocated to the 800 ft. level; the new module included an integrating filter that 
improved the effective energy resolution of the system. 

The start and stop times of each data file run were determined from the 
acquisition computer's internal clock, which is synchronized using the Network 
Time Protocol with a stratum 2 server. These times were used for the primary 
determination of the experiment's live time. In addition, the flash ADC module 
was set to generate periodic sampling triggers every 0.8 s, and these triggers 
were counted for a subset of the data as a cross-check and found to agree within 
0.1%. The time base on the flash ADC board is derived from a crystal oscillator 
that is specified for accuracy and stability to ±50 ppni. 

2. Results 

We have determined that the rate of throughgoing muons on the surface, on 
the lower floor of the SDSTA administration building (~1 m w.e.), is 

(1.149 ±0.017) X 10^2 cm"2 g-i gj.-i ^ 

On the 800 ft. (0.712 km w.e.) level, in the former blasting cap storage area 
near the Ross shaft, the flux is 

(2.67 ±0.06) X 10"^ cm~2 s"^ sr"^ . 



Site/Level 



Surface 



800 ft. 



2000 ft. 



Equivalent depth (km w.e.) 0.001 

Live time (s; ± 0.1%) 3794 



0.712 
6716383 



1.78 
1205647 



Three-fold coincidence analysis: 

Number of events (±^/]V) 15256 

Efficiency of energy cut 0.823 ± 0.009 

Solid angle (sr) 0.457 ± 0.001 



7138 

0.938 ± 0.017 

0.457 ± 0.002 



108 

0.823 ± 0.010 

0.457 ± 0.001 



Flux (cm 



(1.149 ± 0.017) xlO"^ (2.67 ± 0.06) xlO"'' (2.56 ± 0.25) xlO" 



Four-fold coincidence analysis: 

Number of events (±VN) 6930 

Efficiency of energy cut 0.734 ± 0.012 

Solid angle (sr) 0.226 ± 0.001 



3261 

0.837 ± 0.021 

0.226 ± 0.003 



50 

0.734 ± 0.017 

0.226 ± 0.001 



Flux (cm 



(1.184 ± 0.025) xlO"^ (2.76 ± 0.09) xlO"" (2.69 ± 0.39) xlO" 



Table 1: Summary of results from different sites in the mine, showing the consistency of the 
three-fold and four-fold coincidence analysis methods. 



Finally, on the 2000 ft. (1.78 km w.e.) level, at the first transformer pad site in 
the drift between the Ellison and Yates shafts, the flux is 



(2.56 ±0.25) X 10 



"'^ cm ^ s 



Table [T] provides more details on the measurements that led to these flux values, 
and it compares them to the corresponding results from the four-fold coincidence 
analysis method. 

3. Conclusion 

Figure |3] shows the muon flux as a function of the depth, comparing the mea- 
sured fluxes to the prediction of a parameterization model [Ij . The agreement 
between the measurements and predictions for different levels is within 20% in 
all cases. The integrated muon fluxes at the different levels are also compared 
to the prediction of a flat-earth model [4]. Figure [4] illustrates the agreement, 
which is clearly excellent over flve orders of magnitude in flux. We have consid- 
ered the effect of the large open cut in the mine; while it is quite deep, the slant 
depth to our experimental sites is large enough that it has a negligible effect on 
the flux. Consequently, we have used these models to predict the differential 
and integral fluxes at levels of scientific interest, as summarized in Table [2] 
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Figure 3: Muon flux as a function of deptii, placing these measurements in the context of pre- 
vious results from Castagnoli [5], Barrett [6], Miyake [7], WIPP [S], Soudan [9], Kamioka [lO] . 
Boulby [llj . Gran Sasso |12l . Frejus |14| and Sudbury |15| . and comparing them with a pa- 
rameterization model |1] . The three data points for Homestake are from this work (filled black 
circles) and Cherry et al. [3] (filled black square). 



Site 



Flux per unit 
solid angle 



(cm 



Integral flux 



sr ^) (cm -^ s ^) 



4.40 X 10-9 
1.65 X 10"^° 
6.86 X 10-" 



4850 ft. 
7400 ft. 



(1478 m) 
(2255 m) 



3.85 X 10- 
2.21 X 10" 



8000 ft. (2438 m) 1.19 x 10" 



Table 2: Differential (per unit solid angle) and integral fluxes calculated at several levels in 
the Homestake Mine using modeling techniques described in |4| that have been calibrated by 
the results described in this paper. Major experimental campuses are proposed for the 4850 
ft. and 7400 ft. levels, while 8000 ft. is the deepest level. 
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Figure 4: Integrated muon flux as a function of depth, compared against a flat-earth model. 
This model is used to extrapolate to the 4850 ft. and 8000 ft. levels at Homestake. 
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